Transport in Nanotubes: Effect of Remote Impurity Scattering 
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Theory of the remote Coulomb impurity scattering in single-wall carbon nanotubes is developed 
within one-electron approximation. Boltzmann equation is solved within drift-diffusion model to 
obtain the tube conductivity. The conductivity depends on the type of the nanotube bandstructure 
(metal or semiconductor) and on the electron Fermi level. We found exponential dependence of the 
conductivity on the Fermi energy due to the Coulomb scattering rate has a strong dependence on 
the momentum transfer. We calculate intra- and inter-subband scattering rates and present general 
expressions for the conductivity. Numerical results, as well as obtained analytical expressions, show 
that the degenerately doped semiconductor tubes may have very high mobility unless the doping 
level becomes too high and the inter-subband transitions impede the electron transport. 



I. INTRODUCTION 

Carbon nanotubes, discovered last decadei, attract 
attention of physicists not only due to their beauti- 
ful symmetry, extreme quantum scale and quasi-one- 
dimensional (ID) structure, but also due to numerous 
application, already existingSii4 and foreseen in future. 
Among those, the applications, which will use transport 
properties of nanotubes^, are sensitive to details of a 
charge carrier scattering. It was shown that the scatter- 
ing in an ideal metallic single-wall nanotube (SWNT) is 
greatly diminished as in a full symmetry of an armchair 
nanotube the backscattering between states at a Fermi 
level (in a thin shell near an electrochemical potential 
at T 7^ 0) is forbidden^. In contrast, an electron in a 
doped semiconductor SWNT , which has different symme- 
try (zigzag or other chiral symmetry), can be backscat- 
tered. In field effect devices the semiconductor nanotubes 
are used. To obtain high conductivity the doping level is 
controlled, thus, our study of the transport mechanisms 
and the semiconductor SWNT conductivity^ and their 
dependence on the doping level becomes very timely. 

Among most important mechanisms of elastic* scat- 
tering, responsible for the low temperature resistivity 
of nanoscale systems, one may encounter a defect scat- 
tering, a Coulomb center scattering, and an electron- 
electron interaction; the latter is beyond the scope of 
one-electron theory and will be studied elsewhere. First 
two are distinct as a scattering potential is long ranged 
for the Coulomb center, and it is short ranged for a neu- 
tral (mechanical) defect, which may be a substitutional 
(not charged) impurity or another type of a lattice defect. 
The charge carrier scattering by the Coulomb center, 
charged impurity and/or ionized dopant is in the focus of 
our paper. At low temperatures this mechanism is com- 
parable or even stronger than other scattering mecha- 
nisms. Modern high mobility semiconductor devices have 
so small size scale that the device channel may contain no 
substitutional impurity/dopant atom. The remote impu- 
rity scattering may still limit the transport in this case. 



Recent experiments^ showed that even long channel 
semiconductor SWNT devices may have very high mo- 
bility at the high doping level. Our theory addresses this 
case, for the first time taking into account Coulomb scat- 
tering at the remote impurities and including effects of 
an inter-subband scattering at the high (doping) Fermi 
level. 

Specifics of the Coulomb interaction in nanotubes is 
that the effective potential, seen by the delocalized elec- 
trons, has a cut off at the radius of the tube, R. Thus, 
even if the Coulomb impurity is at the closest (atomic) 
distance from the tube surfaceJl, the effective Coulomb 
potential is the potential of a remote scatterer. The situ- 
ation is similar to what is known for HEMT devices. We 
propose the remote impurity scattering mechanism to be 
responsible for a residual resistance of the semiconductor 
nanotube at low temperature. It may define a limit for 
an electron mobility (ON current) in the nanotube field 
effect transistor if other scattering mechanisms are less 
effective. Same theory is applicable to the scattering in 
metallic nanotubes which are proposed for use as inter- 
connects and ballistic wiresi^ and especially in metallic 
field effect transistors 

We develop below a theory of the remote impurity scat- 
tering for the SWNT device, calculate a scattering rate at 
one surface impurity, scattering by a dilute distribution of 
random impurities (in a Born approximation) , a zero and 
finite temperature conductivity in a drift-diffusion (DD) 
model. It is known that in pure (mesoscopic) systems, a 
quantum interference correction to the DD conductivity 
depends on the system dimension and at D = 1 it is of 
the order of the conductivity itself. This manifests as a 
ID localization. However, for a quasi-lD system of the 
nanotube it may not hold due to fast phase breaking. 
It may explain why considerably simple classical one- 
electron theory describes the nanotube experiment fairly 
welli^ii^ii^. Among possible mechanisms responsible for 
the fast phase breaking we notice: e-e and e-phonon in- 
teractions. The latter mechanism was shown to be very 
fast in SWNTs in optical studiesii. We will discuss the 



2 



role of e-e interaction later on. 

We calculate dependence of the SWNT conductivity on 
the Fermi level. Strong degeneracy of carriers in the semi- 
conductor tubes is required to obtain the high mobility 
We show that this is consistent with calculated scattering 
rates. The same condition of high Fermi level allows one 
to apply the Boltzmann equation for transport calcula- 
tions. 

It is commonly used that the conduction in SWNTs 
happens via states of one subband (of orbital quantiza- 
tion). We note that at the high doping level the new 
conduction channel (through the next subband) must 
appear. Then, an interesting question opens: how the 
scattering between different suhhands may change the 
conductivity. In semiconductors role of an intervalley 
scattering is not very important. The difference in the 
transport at high fields consists mainly in the renormal- 
ization of an effective mass of the electron of the con- 
duction band. We will show that for the quasi~lD band- 
structure of the SWNT, where a phase space of scatter- 
ing is very restricted, the opening of a new scattering 
channel may change the conductivity qualitatively, espe- 
cially for the metallic armchair nanotubes and for high 
mobility semiconducting nanotubes. We found that for 
the Coulomb scattering mechanism it results in a giant 
drop of the conductivity at low temperature. We present 
a general expression for the concentration dependence 
of the conductivity and study its analytical limits that 
clearly demonstrates the physics of this effect. 

ID conductivity, ctid, is known to have a dimension 
of a diffusion coefficient, (length)^ /(time) (see for exam- 
ple, Ref.^»). Thus, the conductance of the ID system is 
(Jid/L where L is the system length. We demonstrate 
below that a simple expression is valid for the DD con- 
ductivity of the SWNT: aio ^ GoA, where Go = e^/2TTh 
is the quantum of conductance and A ~ wr is a mean 
free path of the electron {v and r are the electron ve- 
locity and lifetime). We calculated A as a function of 
the Fermi energy and a strength of the random impurity 
potential. Within the model discussed, the conductance 
can be written as: 

here R is the SWNT radius (characteristic transverse 
length of the ID channel), Ep is the electron Fermi level. 
Us is a strength of the ID potential of the Coulomb impu- 
rities, Eg is the SWNT gap, /i and /2 are some power law 
functions, and g{x) is an exponential function of its argu- 
ment. The conductivity (and conductance) is so sensitive 
to the change of the Fermi level (exponential function of 
Ep) because of the exponential dependence of g{x) that 
reflects a strong dependence of the Coulomb matrix el- 
ement on a transferred momentum. Allowed backscat- 
tering transitions within a single subband have a large 
momentum transfer at large Ep^ while the transition be- 
tween subbands may have smaller momentum transfer. 
Thus, an appearing of the new scattering channel dras- 



tically decreases the mean free path of the electron. A, 
which becomes small and G drops several orders of mag- 
nitude. We notice that the power law function fi{x) » 1 
because of its argument is large: Ep ^ Ug. This inequal- 
ity is also a necessary condition for applicability of the 
Boltzmann equation. 

Below we consider charged impurities which are likely 
presented on a surface of any substrate. The deposition 
methods which are currently used in fabrication of the 
nanotube devices may produce such impurities in large 
quantities. Charged impurities are known to exist at the 
surface of Si02 substrate, commonly used for nanotube 
devices. Use of high-e dielectrics may even increase the 
role of this scattering mechanism. 

We assume that the impurities are (single) ionized and 
uniformly distributed on the insulating surface with a 2D 
density Us . Both assumptions are not vital for the model 
and make no qualitative change in the final results. How- 
ever, the derivation of the Eq.(^, for example, is more 
clear in this geometry. The generalization of our theory 
to the case of 3D distribution of the Coulomb centers 
is straightforward. It gives a description for the SWNT 
embedded in an insulating matrix and will be published 
elsewhere^S. For the Coulomb substitutional impurities 
located directly in the lattice of the SWNT and for the 
charged impurities encapsulated inside the tube, an or- 
der of magnitude estimate for the DD conductivity can 
be obtained by substituting [/^ by WgEc, where Ws is 
the strength of a random potential of ID impurities and 
Ec ^ e^/C is a Coulomb charging energy of the tube. 
We notice the logarithmic divergence of the latter with 
the L as well as the same logarithmic divergence of the 
ID Coulomb matrix elemenliSii. 



II. REMOTE IMPURITY SCATTERING RATE 
A. Model assumptions 

The nanotube is situated at the van der Waals distance 
from the surface of the substrate. This distance is about 
3.4 A, and the Coulomb centers are removed from the 
device channel. The Born approximation, which implies 
independent scattering events, is used in what follows to 
calculate the elastic scattering rate. We consider scat- 
tering of the electrons in different subbands and between 
subbands. One can apply this theory to the transport 
in multiwall nanotubes, which are believed to conduct 
by the outermost shell. Here, we restrict ourselves to 
the case of single-wall tubes and consider armchair and 
zigzag SWNTs, though, the final result is more general 
and can be, possibly, used for an arbitrary tube. 

The Coulomb scattering manifests itself in low- 
dimensional systems. The long range Coulomb poten- 
tial is known to be underscreened in a ID caseSitS^, 
in nanowiresi^ and nanotubes2i24. Thus, the Coulomb 
scattering becomes the most important scattering mech- 
anism at certain conditions. It is well known that the 
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random potential results in a localization of carriers in 
an infinite ID systemSSi. We assume that an effective 
phase breaking mechanism exists in the SWNTs, which 
destroys the interference. It is known that the e-e scat- 
tering time is very short for these systems. Even though, 
the e-e interaction in a single ID band cannot suppress 
the localization. Because the total momentum of the 
electronic subsystem is conserved. However, in the nan- 
otubes there is a number of different e-channels which 
may be not coherent. It is clear that the e-e scatter- 
ing between the electrons that belong to the different 
bands breaks the phase of the wave function and, thus, de- 
stroy the interference. What are these different e-bands? 
The transport in nanotubes occurs via the band formed 
by highest valence electrons, so called, 7r-electron band. 
The rest of the valence electrons are localized and form 
low-lying a-bands. The e-e scattering between these 
two different bands does not preserve the electron phase. 
Also, in a real experimental situation, in the SWNT rope 
or in the multiwall nanotube, the scattering between tt- 
electrons at the different walls/tubes destroys the inter- 
ference as well. For what follows we accept that the phase 
breaking time is short enough to neglect the interference 
correction and use the Boltzmann equation approach for 
the calculation of the DD conductivity of the SWNT. 



B. Envelope wave functions 



pseudospin ( distinguishes between states of valence and 
conduction bandSS. 

So far, we considered an equilibrium (non-perturbed) 
electronic wave functions. We assume, as usual, that the 
scattering can be modelled perturbatively if the interfer- 
ence terms are negligible as it was discussed in the last 
section. The perturbation operator is the Coulomb po- 
tential: 



where e < is an electron charge, x, y and z are the 
coordinates of the electron. These three coordinates are 
not independent as the electron motion is restricted to 
the surface of the cylinder, e* is an effective charge of 
an impurity, its position is given in a Cartesian coordi- 
nate system as [Xi, Yi, Zi]. For the scattering at a single 
impurity, a relative position of the Coulomb center along 
the (infinite) nanotube, Zi, may be chosen arbitrary. The 
coordinate Xi (normal to the substrate surface) approxi- 
mately equals h, a negative height of the nanotube, which 
assumes that the impurity size is negligible and that the 
impurity is not buried in the substrate (both assump- 
tions are reasonable but the model works without this 
simplification as well). We will define an effective charge 
of the impurity, e* , in the last section, when discussing 
the screening. 



We use SWNT envelope wave functions^^ for the 
electronic structure calculation, which are obtained as 
a solution of a tight-binding (TB) Hamiltonian for tt 
electronsSi. Our approach is very close to what was pre- 
sented in Refi2& and then widely used in the nanotube 
literature, so we skip details and give only the final wave 
functions for the two-band scheme (tt electrons only): 
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'm,k,c} = ^(1^) + C,c,nk\B))e"- e 



(2) 



here an orbital momentum m labels orbital subbands of 
the SWNT electronic structure, k labels states with a 
longitudinal momentum, both m and k are good quan- 
tum numbers (discrete and continuum, respectively) for 
an ideal, long enough nanotube. C = ±1 is a pseudospin. 
A pseudospinor vector is formed by a two-component 
amplitude of the wave function defined for two atoms in 
a graphite unit cell (A and B). Coordinate along the tube 
is z, and a is an angle along the nanotube circumference 
(by this we explicitly assume the electron to be confined 
to a surface of a cylinder of fixed radius R). The com- 
ponents of the pseudospinor are c-numbers, in general, 
defined up to an arbitrary phase. It is taken such that 
the coefficient for the A-component is purely real and 
equal to 1/ -^2, then a matrix element of the dimension- 
less TB Hamiltonian gives the second component of the 
pseudospinor—. This determines a dependence of c,„fe on 
the subband index, m, and the ID momentum, k. The 



C. Matrix element of the impurity potential 

One needs to know matrix elements of the Coulomb 
potential between the TB wave functions of the elec- 
tron to calculate the scattering. The potential of the 
remote impurity is smooth at the surface of the nanotube 
and, therefore, it is almost constant within the unit cell. 
Hence, the matrix elements of the potential with the en- 
velope wave functions l(2Jl can be approximated by the 
ID Fourier components of Eq.©: 
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Ilm-n\{\k-k'\R)Klra-n\{\k~k'\p) 



here {x) and Kf^ (x) are the modified Bessel functions^S 
(of imaginary argument) of the order /i, r is the vector 
between the impurity center and the point on the surface 
of the nanotube, R is the nanotube radius. In a cylin- 
drical coordinate system, ip is the angle of the impurity 
position, p is the distance from the axis of the nanotube 
to the impurity and Z is its longitudinal coordinate. 

These matrix elements of the potential are needed for 
calculating the remote scattering rates. Besides that, the 
analytical expression for the remote potential is in- 
teresting by itself. We are not aware of a calculation 
of the Coulomb potential for the charge center removed 
from the nanotube. This formula gives a generalization 
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of an expression for the ID Fourier transformation of the 
Coulomb interaction between charges which are both on 
the nanotube surface (which may be found, for example, 
in Ref.'^^). Let us present analytical limits of our result at 
large and small transferred momentum and demonstrate 
how the dimension of the nanotube system shows up. 

The interaction strength decreases rapidly with the 
transferred momentum q = \k — k'\, which is well known 
property of the Coulomb potential at any dimensions. 
Using the asymptote of the Bessel function at g ^ 



P 



,R \ we reduce the Eq.Q) to 
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(5) 



and recover formula e^^^/g, the expression for the 
Fourier component of the Coulomb potential in 2D^^. 
This is not surprising, because in the short-wavelength 
limit one restores the planar geometry when the curva- 
ture of the graphite sheet becomes unimportant. 

On the other hand, at small q and m — n the ma- 
trix element Q logarithmically diverges as ^ \og{qp/2), 
which is according to ID electrostatics^i. In contrast, 
at m 7^ n and small q, the limit of the matrix element 
has no dependence on q. Instead it is proportional to 
I, and decays exponentially with the trans- 
ferred angular momentum |m — m'| as in the multipole 
expansion series. It is consistent with understanding of 
as a minimum cut off momentum in the Coulomb 
matrix element. 



not coherent between events. Thus, one has to sum prob- 
abilities of the scattering over the realization of impu- 
rities (to be averaged later). Let us write the partial 
probability of the single scattering event using the Fermi 
golden rule: 



9 ^ ^ (6) 



where q is the transferred momentum, and for the SWNT 
the sum has only several terms (less than four) if any, 
which depends on m, n and k quantum numbers. Here, 
h and Y are the Cartesian coordinates of the impurity 
(without loss of generality, we chose the coordinate ori- 
gin such that Z = 0). We remind that the axis of the 
nanotube is at the distance X — h from the substrate. L 
is the tube length. 

The Coulomb centers are distributed on the surface of 
the substrate randomly. One has to perform averaging in 
the plane to obtain a statistical description of the scat- 
tering. We assume that the impurity positions are not 
correlated. Then for the electron with the momentum k 
in the subband |m), an elastic lifetime due to the remote 
impurity scattering is written as follows: 



r """(to. A:) 



dEr, 
dk 



k=qk 



Gjqk) 



(7) 



D. Averaging of the random potential 

In the Born approximation each scattering event is sta- 
tistically independent and the electron wave function is 



where Hs is a 2D density of the surface impurities and qk 
are the solutions of the equation E^^k+q — ^-m.fe- As we 
noticed before this equation may have up to 4 solutions 
within the first Brillouin zone (see Fig^. For example, 
for the armchair nanotube we have: 



qk = — fc ± 2arccos 
and for the zigzag SWNT: 
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We use the notation G{qk) — I^m-n\ (QkR) ^\m-n\ {qkh) for a reduced Coulomb matrix element in Eq.lO. Here the 
factor J- comes after averaging in the plane and equals 



dtK: 
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dtK2n{t) 
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where L are the modified Struve functions-- (of imaginary argument 
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Now we analyze these expressions in the limit of small and large g^. For small qkh one can write: 



Qk ^S{qk) = Qk ^lf,n^n\i<lkR)^[m-n\iqkh) 



TT^/iqk 

ttRR 
16 h 
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if n = m 
if \n — ?7i| = 1 

n — m\ > 1 



r 



(10) 



The numerical factor is rather small in case of inter- 
subband transitions n ^ m and decreases with the sep- 
aration between subbands, |m — n|, rapidly. Thus, the 
scattering into the same subband (transition with the or- 
bital quantum number conservation) is the most effective 
scattering channel for small qk- 



For large qkR the scattering rate exponentially de- 
creases with the transferred momentum ~ exp[— 2(2/i — 
i?)(7fc], due to the exponential decay of the modified Bessel 
function. The main term of a Poisson series of the matrix 
element does not depend on |to — n|: 
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(11) 
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We will use these analytical expressions (|1(JII and for 
the calculation of the scattering rates in what follows. 



III. CONDUCTIVITY: DRIFT-DIFFUSION 
APPROACH 

In this section we calculate the conductivity in a drift- 
diffusion (DD) model which is widely used for description 
of the transport in multiwall nanotubes. It is also appli- 
cable for very long SWNTs if the phase breaking and/or 
inelastic scattering lifetime is short as discussed in Intro- 
duction. 

It is known, that a Schottky barrier forms near the 
metal contacfei^iS. For the short channel SWNT device 
this Schottky barrier almost determines important trans- 
port characteristics. The theory of the nanotube trans- 
port, taking into account phenomena in the contact re- 
gions, is presented in Refs.^^''^^ and we do not address 
this subject in our paper. Instead we focus on the scat- 
tering in the rest of the tube. The contact region has a 
finite length which is about a typical screening length, 
e.g., a distance to the backgat o'^^i^^ . If this distance is 
much smaller than the length of the nanotube, one may 
define the device channel conductivitji^LSI. For the sake 
of clarity, we restrict ourselves to the case of the armchair 
or zigzag SWNT. The generalization of our model to the 
case of any chiral SWNT is straightforward. 

The conductivity in a single channel is as follows: 
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2tt 



e^T{E) dE df 
h? dh 'dE 



dE. 



(12) 



Here g is the degeneracy of the current channel. By the 
channel of conductivity we understand here any fixed 



subband of the orbital quantization which can carry the 
current. t{E) is the transport lifetime; f{E) is an equi- 
librium distribution function. The derivative of the dis- 
tribution function, df /dE is peaked at the electrochem- 
ical potential (delta-function of the Fermi level E = Ep 
for T — 0). However, we keep the integral sign even 
at T = because of several channels corresponding to 
several non-zero terms in the conductivity, as given by 
Eas. H8al8b() . The a has non-trivial temperature depen- 
dence due to the strong dependence of the lifetime on the 
electron energy. 

It is important to include the inter-subband scattering 
terms (if corresponding transition is allowed) because the 
Coulomb scattering rate decreases with the transferred 
momentum (see also Fig. O. With increasing Ep the 
scattering into the same subband may become less ef- 
fective than the scattering into the other subband. We 
demonstrate below that this is the case for the semicon- 
ductor SWNT at the high doping levels. 

The Eq.jTll can be conveniently rewritten as 



^{m, k) 



ml 



,qk ^Q\m-n\ {qk) , (13) 



with use of notations: Us — (ee*) y^nj for a characteristic 
energy of the Coulomb disorder, and Vm = h~^dE/dk for 
the electron velocity. The prime sign reminds that the 
summation is over the roots of Eas. (|8al8b|l . 



Let us first consider the scattering of the electron in 
the same subband \m,k) — > |m, /c -f g). 
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FIG. 1; Electron energy dispersion for armchair [10,10] and 
zigzag [17,0] SWNTs. Inset: Zoom out of the lowest subbands 
and the electrochemical potential (green dashed line). 



A. Armchair nanotubes: Intra— subband scattering 

The expression for the scattering roots for the armchair 
SWNT, Ea. (|8al) . in the limit of intra-subband scattering 
(m = n) reduces to: 



= —kzL 2arccos 



TOTT k 

cos ^ + cos - 



-2fc, 0, (14) 



The last root means no scattering because all quantum 
numbers are conserved and has to be discarded. The first 
root is, in fact, the inter-subband backscattering near 
the same Fermi point (see Fig Q Right Inset), which is 
forbidden for the two crossing subbands of the armchair 
SWNT by the pseudospin conservation ruleS. The root 
q — —2k is the backscattering of the electron within the 
same subband to other Fermi point (see FigQ Right). 
Only this transition between two Fermi points with a 
large transferred momentum {qk ~ in/SV^b for small 
Ep) is allowed by the symmetry of the armchair SWNT. 
Thus, the analytical limit can be used for calculating 
t{E). The transport lifetime of the electron contains a 

"47r(2/i-fl) 



large exponential factor oc exp and, hence, 

this mechanism gives the negligible scattering rate. 

Because of the exponential dependence of the matrix 
element on the transferred momentum, the conductivity 
of the armchair nanotube increases dramatically, when 
the electrochemical potential reaches the next subband. 
This subband has a different symmetry^^ and the inter- 
subband scattering near the same Fermi point (see Fig^ 
Right Inset) is not forbidden. Then, the inter-subband 
scattering rate of the armchair nanotube is high and 
given by an expression which is similar to the case of 
zigzag nanotube. We will consider these two mechanisms 
together in Sec llll Cl 



B. Zigzag nanotubes: Intra-subband scattering 

In the last section we show that the remote impurity 
scattering is negligible in the armchair SWNTs. The 
scattering due to any other long range potential, except 
for torsional phonon modea^S, was shown to be ineffec- 
tive as well^^. The case of the zigzag nanotube is less 



trivial. We substitute Eq.(|7I) into Ea. ((T^ and obtain the 
conductivity at zero temperature 
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Q\m-n\ i\kFm — kpn) 
Vm'^n\kFm ~ kpn] 



(15) 

here g — A is the current channel degeneracy (for the 
spin and orbital momentum); Vm is the electron velocity 
at the Fermi level (near the F-point), Vm = hr^dE / dk. 
Taking into account the non-parabolicity of the energy 
dispersion of the to*'^ subband, -Efc.m we obtain: 



667 
AhEp 



EI-EI^Vf, 



'1- 



E^ 

m 

E^ ' 



(16) 



here Vf = 3bj/{2h) is the Fermi velocity in the (metallic) 
massless subband. The Fermi momentum of the electron 
in the m"^ subband is: 

kmiEp) 



2 
367 



El 



E^ 



(17) 



where the bottom of the m*'' subband is 



Em = 7 



1 + 2 cos 



N 



hVprn 
3i? ' 



(18) 



These expressions are essentially similar but not equiva- 
lent to ones obtained with k ■ p method because of non- 
parabolicity of the energy dispersion in tight-binding 
model. 

For the electrochemical potential located within the 
lowest subband, Em < Ep < Em±i, the only one level 
of orbital quantization is populated at T = 0, which has 
m= {N + Mod3[N])/3. 

Because the inter-subband transitions become allowed 
only at Ep > Em±i the single backscattering term with 
m — n has to be substituted in the Ea. (|15l) which gives 



^J-^ = GogVmT{m,kp) 



Gggh^ 
8C/2 



2kFv'ig-\2kp) (19) 



here Go — 12'Kh is the conductance quantum, and the 
electron velocity is 



"El 



if ?7 ^ Mod3(iV) 7^ 
if ry = Mod3(iV) = 



(20) 



here the index 77 = Mod3(A^) distinguishes between the 
zigzag metallic (77 = Mod3(A^) — 0) and semiconductor 
SWNTs (?7 = Mod3(iV) ^ 0). 

Let us now apply the expression H1U|) for analysis of 
the conductivity at small Ep. Because at small kp the 
function Q^^{2kp) does not depend on kp in the leading 
term, we find the DD conductivity of the zigzag tube at 
the low doping level depends on the concentration as: 



.(1) 



GogEphVp ^ 
4C/2 ■ 



3/2 



if 
if 



77 = 0. 



(21) 
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We drop the argument of in the expression above be- 
cause its hmit is A/tt'^ at small kp. The DD conductivity 
is different for the metallic and semiconductor tubes: it 
is linear in Ep for the metallic (zigzag, 77 = 0) SWNT 
and it depends on the Fermi energy as 

for the degenerately doped semiconductor SWNT 
(zigzag, 11^0). 

The conductivity at small Ep depends linearly on Ep 
owing to the matrix element of the Coulomb potential of 



the 2D remote centers. The square of the matrix element 
is ~ , because in our model the centers are distributed 
on the (2D) surface of the substrate (the generalization 
to the 3D case is obvious). The linear dependence of the 
square of the matrix element in results in the linear 
dependence of the conductivity on Ep. 

At larger Ep the energy dependence of the transport 
lifetime (of the matrix element of the transition) is dif- 
ferent. In this case, the Ea. Hll(l has to be used, which 
results in a fast exponential growth of the conductivity 
because of the large suppression of the transitions with 
increasing momentum transfer. Then, the Ea. (|21(l trans- 
forms into: 



r(l) 



SRGogEl 

uT 



hEp 



£2 



\ TThVp 



^ El V El 



Ei 



exp 



A{2h-R) 



Ep 



£■2 



hVp V ^ El 



here ep — Ep/Em is the dimensionless Fermi level of 
the semiconductor tube: tp ^ \ and ep < 1 in the ex- 
pression (|22|l . This equation gives the conductivity of 
the zizgaz nanotube for the transport through only one 
channel (of the lowest subband). The result of our calcu- 
lation (the Eqs.(j2U and (|22|l l is presented in the Figure 
121 where the drift-diffusion conductivity of the zigzag 
SWNT is plotted as a function of the Fermi level (dop- 
ing level). With increasing doping level the conductivity 
grows exponentially. This is not because of more carri- 
ers are available for transport but due to the energy de- 
pendent transport lifetime that grows exponentially with 
the increase of the momentum transfer between initial 
and final states. This momentum transfer is 2kp oc Ep 
when the linearization of the TB Hamiltonian is possi- 
ble. The momentum transfer increases until the Fermi 
level reaches the next subband. Then, suddenly, the new 
backscattering channel opens. The momentum trans- 
fer between the next subband and the lowest subband 
is small at this critical doping level. Thus, the remote 
impurity scattering becomes very efficient and the con- 
ductivity drops several orders of magnitude. 

The maximum conductivity may be reached when the 
Fermi level is close to the edge of the second subband 
but lower than it. We note that the phonon scattering 
may limit the transport, at least, at high temperatures^. 
Therefore, for not very low temperature and high enough 
doping level, it may be possible to observe switching from 
the remote scattering to the phonon scattering mecha- 
nism. 

Also at r 7^ 0, the finite temperature distribution func- 
tion makes the a{Ep) dependence smooth near the step 



(22) 

l)''^exp [%{2h~R){ep~l)] 
I 




0.4 0.6 0.8 1.0 1.2 



FIG. 2: Conductivity of a zigzag [17,0] SWNT vs. the elec- 
trochemical potential (the doping level). 



at ei? '--^ 1, as it will be shown next. 



C. Inter— subband scattering 

The remote impurity scattering between the subbands 
may happen only when the doping level is high enough 
to essentially populate the second subband, Ep > Em±i- 
Then, the scattering rate becomes high and the mean free 
path becomes short. The expressions for the DD con- 
ductivity of armchair metallic, zigzag quasi-metallic and 
zigzag semiconductor nanotubes are essentially the same 
in this region. This is because the DoS in the vicinity 
of the Fermi point is a universal function^ of the energy 
(doping level), and because there is no special selection 
rule for the transitions between different subbands of the 
orbital quantization. Let us consider the single Fermi 
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point and find which scattering channel defines the con- 
ductivity at Ep > Em±i. As before, because of the large 
momentum transfer, we neglect transitions between dif- 
ferent Fermi points, which are possible for the armchair 
SWNT. As shown in Fig|31 we have two left-going and 
two right-going (current) states/channels (to be multi- 
plied with the spin and orbital momentum degeneracy 
g = 4). We introduce four scattering rates Tj"^ and cal- 
culate it with the Ea. ((n|) . 

Let us assign the index z = 1 to the subband with 
the largest momentum at the Fermi level, kpt (which 
may change with changing the Fermi level if the sub- 
bands cross). The intra-subband lifetimes are 
given by terms of Eq.Q with n — m. There are two 
inter-subband lifetimes: r^^^' and r™^' for the transi- 
tion with the same/opposite sign of the electron velocity 
in the initial and final states \n > and |m |n ± 1 >. 

For certain tube symmetry the subband crossing may 
happen with increasing Ep. When the Fermi energy is 
above the crossing point the subbands 1 and 2 are in- 



and we study below the limiting cases where simpler an- 
alytical expressions are available. 

At small kp2, at the second subband edge (see Fig|31 
the Fermi level is at lower/pink line), the momentum 
transfer of the inter-subband scattering (1 ^ 2) is larger 
than of the intra-subband scattering (2 ^ 2) (see Figure 
13 Inset) and the last scattering channel is more effec- 
tive. The corresponding contribution to the cr is oc W2T22, 
where the velocity at the Fermi level is given by Eg. lit) |) 
(the lowest subband i = 1 has the larger velocity). 

According to inequality: 

> r-/ > > t',' (24) 

the most important term in the cr is cx viTr, if the shortest 
time is T22- Thus the DD conductivity reads as 

a = GoQViTr -\ ~ g°^2 ^1^2|fcl - ^2! Q~^{\kl - fc2|), 

' (25) 

where vi = Vp^l - {hVpf / [iREpY and V2 = 
Vp^l - {2hVpY/{ZREFY. 

With increasing Ep (see Fig|31 Fermi level at up- 
per/green line), since (722 > '?12, then, r22 < tv (see Inset 



terchanged in the equations given below. The case of 
the metallic (armchair or zigzag) SWNT is similar to 
the case of noncrossing subbands and will not be con- 
sidered separately. There is only one difference for the 
armchair SWNT as compared with the zigzag SWNT 
case: the dispersion in the subband 1 is massless. Thus, 
vi = Vp is the constant and, by symmetry, rn = 00 
(" intra-subband" scattering in the lowest subband is not 
allowed) . 

We recalculate the distribution functions in all 4 chan- 
nels using Boltzmann equation and taking into account 
the inter-channel transitions. When the scattering is 
weak one can neglect it and use the equilibrium distri- 
bution function. This is not the case for the SWNT at 
ep > I, where the scattering rates at the edge of the 
second subband are very high. 

General expression for the DD conductivity with the 
non-zero inter-subband scattering is rather cumbersome 
even in the approximation of two closest subbands: 



(23) 




FIG. 3: Bandstructures for zigzag SWNTs of two types. 
Left: [3q-hl,0] SWNT; Right: [3q-l,0] SWNT. Doping level 
is shown as horizontal lines; pink/green (lower/upper) line is 
for low/high doping level. Insets show the scattering rates for 
different scattering mechanisms as a function of energy as in 
left and right diagrams respectively. 



of Figure OJ, and the leading term of Ea. (|23|l is due to 
the intra-subband transition (2 2). Then, for the arm- 
chair and zigzag (77 7^ -1-1) nanotubes, the main term in 



a = Gog- 



, {Vi+V2f 
VlTr H r22 

2vi 



(Wl - V2f Tr ^ Tr 

n T22 + V2 T22 

2^1 Tm Til 



2V2 T22 , 2Tr , V2 T22 



1 



Trn Vl T„i Tn Vl Tu \ T, 



8C/2 



[viV2\ki - k2\g-H\ki - k2\) {1- 



\ki-k2\g-'{\ki-k2\ 
\ki + k2\g~^i\ki + k2\ 



(t;i+^^ 1 / , {vi-v2r \ \ki - k2\g~H\ki - k2\) 
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the DD conductivity is 



(T, cm/kOhm 



^ {Vl + V2? 



(26) 



8C/2 



2vi 



2k2 g~\2k2) 



In the zigzag 77 = +1 SWNT, the crossing of the lowest 
subbands occurs in the studied region (shown in Fig. O 
Left). In that case kpi > kp2 for the Ep is lower than the 
crossing point, and kpi < kpi for the Ep is higher than 
the crossing point. So the greatest scattering rate is due 
to (2 ^ 2) intra-subband transition below the crossing 
point and (1 1) intra-subband transition above the 
crossing point (see Left Inset of Fig. |2J). Thus, the indices 
1 and 2 must be exchanged in the Ea. (|23|) and Ea. (|26|) . 

By considering these two cases we cover all possible 
situations and present possible analytical expressions for 
the DD conductivity within the remote center scattering 
model. 



IV. CONDUCTIVITY AT THE FINITE 
TEMPERATURE AND IMPURITY POTENTIAL 
FLUCTUATION 

In the last section we studied the conductivity of the 
nanotube in the zero temperature limit. The tempera- 
ture dependence adds to the above result via substituting 
Fermi-Dirac distribution functions in Ea. (|12|l instead of 
step functions as we implicitly used before. 

We present the numerical result on the temperature 
dependence of the conductivity in Fig^ 

We already discussed an important assumption of our 
model: the phase breaking time has to be short enough 
which seems to fulfill for the nanotubes due to the e-e 
interaction and scattering of the electron into the chan- 
nels/bands which are different from the transport chan- 
nel/band. This allows one to neglect the interference 
correction which is normally dominating in pure ID sys- 
tems. However, the bare ID Coulomb potential may still 
localize the carriers in a nanowire for the infinite channel 
length. We remind that the quasi-one-dimensional po- 
tential created at the nanotube channel by the random 
distribution of remote impurities on the substrate surface 
has the cut-off length and, thus, a maximum amplitude 
of the potential. The localization length in a finite sys- 
tem is defined by the average fluctuation of the random 
potential. Our theory is applicable only in the limit of 
the electrochemical potential which is much higher than 
the average fluctuation of the random potential of the 
remote scatterers. 

Since the Coulomb centers are located on the substrate 
and the ID electron is on the nanotube, there are two dif- 
ferent types of averaging for the 2D distribution of impu- 
rities and for the ID random potential for the electron. 
The operator of the Coulomb potential is given by the 




(J, cm/ Ohm 




0.3 



0. i 



0.5 



0.6 0.7 
|Li,eV 



FIG. 4: Temperature dependence of the conductivity of a 
zigzag [17,0] SWNT vs. the electrochemical potential in a 
vicinity of the second subband edge. 



Eq.Q. It creates the ID potential along the nanotube, 
which reads as: 



UiD{z) = Y.V{Y„Zi 



(27) 



where the sum is over impurities that have random posi- 
tions. 

To calculate the average fluctuation of this potential 
along the nanotube we average it over the impurity po- 
sitions: 



= \J {UId) - {UidY = 2Y^2 + 21og2 + log22 ee*V^. 

(28) 

At low temperature, for random impurity distribution, 
this average fluctuation gives an estimate for a lower 
bound of the electrochemical potential at which the Born 
approximation for the scattering is applicable: 



Ep > Us 



(29) 



It is known that the condition of applicability of the 
Boltzmann equation is that the Fermi energy must be 
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much larger than the inverse scattering time. If we esti- 
mate the latter as H/t ^ 411'^ /Ep then we arrive to the 
condition which is similar to Ea. H29|l . 



V. EFFECTIVE IMPURITY CHARGE AND 
SWNT DEPOLARIZATION 

Though, the actual charge of impurity is not known, 
we assume it to be an elementary charge e. However, the 
substrate polarization results in a reduction of this value. 
As long as the distance between the charge center and the 
substrate surface is much smaller than all other lengths 
of the problem: R, h, etc., one may use an effective di- 
electric function of the substrate to define the effective 
charge as e* ^ 2e/(e -I- 1) where e is the dielectric func- 
tion of the substrate and unity stays for the permittivity 
of the vacuum. It is the leading term of expansion series 
of the image charge potential, which has to be kept in 
the remote scattering calculation. 

One must take into consideration an effect of depo- 
larization of the Coulomb potential due to the screening 
by carriers in the nanotube. This changes the remote 
scattering potential essentially. We calculate this effect 
using a continuum model for SWNT electrostaticsS^*^. 
Within the model, the depolarization of the SWNT at 
the distance Dg from a conducting gate is given by the 
following expression: 



1 



eff 



-V = 



c; 



(30) 



where the depolarization a is written in terms of the 
geometric capacitance per unit length of the SWNT 
= 21og(2Dg/i?) and the quantum capacitance per 

unit length, which equals Cq^ = l/^e^VM) for the metal- 
lic and degenerately doped semiconductor tube. 

The depolarization can be taken into account together 
with the substrate image charge effect simultaneously: 



1 Cn' + Cr 



(31) 



VI. CONCLUSION 

In summary, we have developed a microscopic quan- 
tum mechanical model of an electron scattering by re- 
mote Coulomb impurities lying on the substrate surface. 
We consider a uniform 2D distribution of the charged 
impurities. Theory is readily generalized for the three- 
dimensional distribution of the impurities which case is 



applicable to the modeling of an insulator surrounding 
the nanotube channel (to be published elsewhere). A gen- 
eral expression for a Fourier transform of the Coulomb 
potential of a charge which is removed from the nan- 
otube surface is given. The matrix element for calculat- 
ing the remote impurity scattering in the nanotubcs is 
derived. We obtain a zero and finite temperature con- 
ductivity within a drift -diffusion transport model. We 
assume that the e-ph and e-e scattering between the elec- 
trons in the 7r-electron band (transport band) and the 
electrons in other bands (or/and other tubes in the rope, 
or/and other walls in a multiwall nanotube) presents a 
fast mechanism of the phase breaking. It allows us to use 
the Boltzmann equation and neglect interference correc- 
tions for the classical conductivity. Selfconsistent calcu- 
lation of the SWNT depolarization factor, taking into 
account image charges in the substrate, is performed and 
yields an effective charge of the impurity, used for com- 
puting the scattering rate. 

The remote impurity scattering is almost negligible for 
the armchair SWNTs if the Fermi level is below the edge 
of the second subband, which is consistent with other 
calculations for other scattering mechanisms with long 
range potentials. In contrast, the scattering rate for the 
zigzag SWNT is high enough. In general, the DD con- 
ductivity of the SWNT is ~ GqA, where the mean free 
path, A ~ VT. Analytical expressions for the conductivity 
are obtained in the limit of small and large momentum 
transfer. We studied numerically dependence of the con- 
ductivity on the Fermi level position and found that the 
highest conductivity (of a semiconductor SWNT) may 
be observed for the degenerately doped SWNT when the 
Fermi level is close to (but lower than) the second sub- 
band edge. The dependence a{Ep) is not monotonic. At 
even higher doping level, the conductivity is low when we 
take into account the inter-subband scattering. This is 
because A becomes short with opening a new scattering 
channel, when the Fermi level is higher than the second 
subband edge. Thus, over-doping of a semiconducting 
nanotube does not improve its transport properties. Our 
theory may be applicable for multiwall nanotubes, al- 
though, an additional analysis will be required, especially 
because of different screening. 
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